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AB CT 


This thesis describes the steps necessary to develop an acoustic vorticity meter for the 
atmosphere. The analysis is based on Benthic Acoustic Stress Sensor (BASS) technology 
that is currently used for similar acoustic measurements in the ocean. Compared to sonic 
anemometer measurements, the BASS measurements of velocity are not only made in a 
different fluid but in a different way. Due to these differences, the physical make up of 
BASS needed to be altered, and the validity of the measurement technique had to be 
explored. 

The alterations to the BASS hardware occurred for several reasons. Because 
attenuation of sound is much higher in air than in water for the same frequencies, it was 
necessary to change the transducers. The generally faster and unidirectional mean flows 
that are present in the air encourage open measurement volumes which the BASS vorticity 
meters do not have. The difference in group speed of sound is different for water and air, 
and this forced a change to the timing and burst generation board of the BASS vorticity 
meter. 

The measurement technique used by the BASS instrumentation is validated by the error 
analysis in the text. Because the BASS instrumentation actually provides a time 
difference, the equation used by the BASS instrumentation to compute velocity was 
assumed throughout the error analysis. The error analysis shows that the combination of 
BASS measurement techniques with a temperature sensor will provide errors that are less 
than 2% of the velocity. 

The types of measurements that an atmospheric vorticity meter would provide to a 
researcher are described in the text to show the meter’s potential. If deployed on a buoy, 
a vorticity meter could measure shearing of the wind close to the surface of the waves. If 
deployed at heights much greater than its path lengths, an atmospheric vorticity meter 
could provide three-dimensional vorticity measurements which would provide a unique 
measurement of a fundamental characteristic of turbulent flows. 
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СНАРТЕК 1 


INTRODUCTION 

The turbulent motions in the oceanic and atmospheric boundary layers are primarily 
responsible for the transport of momentum, heat, gases, and mass toward or away from 
the interface. If one wishes to study the turbulent processes that occur on either side of 
these interfaces, fast response instrumentation is required. Over the past several decades, 
sonic anemometry has advanced to the point that sonic anemometers are routinely used to 
measure turbulent fluctuations in the atmosphere. Several types of sonic anemometers 
determine the instantaneous wind velocity by measuring the time of flight of acoustic 
pulses in air (Schotanus et al, 1983). Sonic anemometers record the times of flight of 
separate acoustic pulses in air in opposite directions along each acoustic path. These 
measurements are made consecutively. The high frequency of measurement allows the use 
of sonic anemometers for measurements of turbulent motions in the atmosphere. When 
these measurements are taken near the air-sea interface, we can gain insight into the 
exchange between the two fluids. 

The BASS, Benthic Acoustic Stress Sensor, velocity meter is frequently used to fulfill 
a similar role in the oceans to that of sonic anemometers in the air. The BASS was 
originally designed to operate in water near a fluid-solid interface while sonic 
anemometers operate in air near a fluid-fluid or fluid-solid interface. The BASS 
instrumentation and structure have recently been modified to create a family of three axis 
acoustic vorticity meters for use in water. These vorticity meters have been deployed to 
measure the shear in the upper ocean boundary layer, internal waves in the thermocline, 
and turbulence in the bottom boundary layer (Thwaites et al., 1995). 

The measurement of vorticity as a variable would greatly assist atmospheric turbulence 
research. Therefore, a group of air-sea interaction researchers at Woods Hole 
Oceanographic Institution (WHOI) decided that it would be extremely useful to transfer 
this acoustic vorticity technology to the atmosphere. An atmospheric vorticity meter 
could be used to make vorticity measurements at the bottom of the atmosphere that are 


similar in nature to those made by its aquatic counterpart. It could also be used to 





measure the shear near the surface over ocean waves or over land. 

The BASS measurements are different from the sonic measurements in several 
important ways. The BASS uses transducers that are designed to make measurements in 
the ocean and not ın the atmosphere. The BASS timing circuitry produces acoustic pulses 
that travel in opposite directions along an acoustic path and occur concurrently rather than 
consecutively. Therefore, the BASS measurements along each path are time differences 
rather than two recorded times of flight per path. Due to these differences in the 
measurement technique, the transfer of a three axis acoustic vorticity measurement to the 
atmosphere will require some modifications and additions. 

In order to operate in the marine boundary layer, an atmospheric vorticity meter must 
be rugged and reliable enough to withstand the harsh conditions encountered over the sea. 
Another objective in the development of the atmospheric vorticity meter is to give it a 
large dynamic range. With a large dynamic range, it will be possible to use the vorticity 
meter in a wide spectrum of conditions, which would lead to broader use of the 
instrument. It must also be portable with path lengths of about 0.15 meters in length. 
With these built in characteristics, the atmospheric vorticity meter should be adequate for 
shipboard or buoy use. At sea use is the main target for the final stages of development 
because less is known about boundary layers over the oceans than over land. This is 
particularly true in the layer very close to the sea surface where wave induced motions 
complicate the flow field. 

The motivation for measuring vorticity is that in many ways it is the superior estimate 
of the turbulent kinetic energy in comparison to an anemometer’s covariance 
measurements. For example, on days with strong winds, shear is the largest contributor to 
turbulent kinetic energy near interfaces. An acoustic vorticity meter is able to give a direct 
measurement of the instantaneous shear. An acoustic anemometer is not capable of 
providing this measurement. 

The sections which follow contain a description of the efforts to modify the BASS 
vorticity meters for atmospheric measurements. Chapter 2 describes aspects of vorticity 
and gives examples of vorticity’s relevance in turbulent studies. Chapter 3 describes the 


physical modifications to which the BASS instrumentation was subjected during its 





conversion for atmospheric applications. Chapter 4 describes some of the errors that will 
be inherent in the initial atmospheric design. Chapter 5 describes the calibrated windtunnel 
that will be used to determine the accuracy of BASS instrumentation in air. Chapter 6 
contains the conclusion. In the next chapter, I begin my discussion with a description of 


the quantity being measured, 1.e., vorticity. 
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CHAPTER 2 


ATMOSPHERIC VORTICITY AND CIRCULATION 

Atmospheric vorticity can be described in several ways, but it can be visualized as a 
vector field which provides microscopic measurements of rotation at every point in the 
atmosphere (Holton, 1979). A vortex can be thought of as a three-dimensional eddy with 
a principal axis. These eddies are the irregular, or turbulent, motion responsible for the 
turbulent transport of momentum, heat, gases, and mass. As an example of the generation 
of eddies, a shearing of the mean wind will tend to produce a vortex with a principal axis 
that is aligned with the mean strain rate (Tennekes and Lumley, 1980). These eddies 
produce effects that are visible to the human eye, such as catspaws on the water. The 
wavelike motions in wheat fields also mark the passing of eddies. These examples 
describe eddies that are associated with forced convection. Forced convection is a 
situation in which shear production of turbulent kinetic energy exceeds buoyant 
production. Forced convection typically dominates on windy days with little surface 
heating. | 

The scintillation that is visible over deserts 1s another example of eddy generation. 
However, this form of eddy generation is associated with free convection. Free 
convection means that buoyant production of turbulent kinetic energy exceeds the shear 
production. For example, free convection occurs when there are light winds with strong 
surface heating. Free convection results in the production of buoyant plumes of air and 
associated vertical, turbulent motion. 

Formally, vorticity is defined as the curl of the velocity vector, which is a measure of 
angular rotation about a local position. This leads to the concept of length scales that 
affect the local position. The maximum length scales that effect the vorticity relate to the 
largest eddies that are important to the local position. For example, in the atmospheric 
surface layer the largest eddies logically scale with the height above the surface layer. 

Mathematically, the local position collapses to a point as the spatial interval lengths of 
velocity measurements approach zero; 1.e., the formal definition of vorticity relates it to 


the microscopic rotation about a point in space. In practice, velocity measurements are 
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separated by some distance, and an average vorticity is measured. This area averaged 
vorticity is closely related to the concept of circulation, which is the measure of the 
macroscopic circulation over some area (Holton, 1979). The circulation over a 


macroscopic path is represented by Stokes Theorem: 
C= $Vedl=[ [(VxV)-dA (2-1) 


where C is the calculated circulation, V is the measured velocity, di is the acoustic path, 
and dA is the area enclosed by the acoustic path. This circulation information is used to 
estimate the area average vorticity over the macroscopic area by using the following 


relationships: 


me 
Sn pev (2-2a) 


(2-2b) 


n 
> | С 


where & is the vorticity, & is the area average vorticity, and A is the area enclosed by the 
sonic paths of the vorticity meter. 

Since time-averaged, area-averaged, horizontal vorticity equals time-averaged, area- 
averaged, vertical shear (Thwaites, 1995), it is useful to determine an expected vertical 
shear to determine the accuracy requirements of a vorticity meter. If a neutral boundary 


layer is assumed, the shearing of the mean wind can be modeled by: 


==. (2-3) 


where u, is the friction velocity, k is von Karman's constant, and z is the height above 


[2 





the ground. If a reasonable friction velocity of 0.20 m/s were assumed at a measurement 
height of 1.00 m, a 0.15 m path length would require a velocity measurement error that is 
at least an order of magnitude less than the 7.5 cm/s velocity shear for the measurements 
to be useful. If this condition could be met, it would be possible to recover the shearing of 


the wind in the vertical from the vorticity measurements. 
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СНАРТЕК 3 


MODIFICATION OF BASS 

The Benthic Acoustic Stress Sensor, BASS, was designed for use in the water, such 
that the BASS instrumentation is not precisely compatible with use in the atmosphere. 
Basic modifications to transducer design, timing and burst generation, and physical shape 
of the sensor are necessary to convert a BASS vorticity sensor to atmospheric use. These 
modifications are necessary because of the difference in how air responds to high 
frequency sound waves. Additionally, the effects of flow distortion by the vorticity meter 


are expected to be different in air versus water. 


INSTRUMENTATION 

A common instrument used for atmospheric boundary layer research is the acoustic 
anemometer. A singular velocity measurement along one acoustic path is similar in nature 
for an acoustic anemometer or vorticity meter, such that an acoustic anemometer suffers 
from the same problems associated with high signal loss through air as a vorticity meter. 
Therefore, it was useful to look at an anemometer during the development of our vorticity 
meter. The advantages that are found in the anemometer should be incorporated into the 
design of the final atmospheric vorticity meter. 

It is common for acoustic anemometers to use open measurement volumes, fifteen to 
thirty centimeter path lengths, and piezoceramic transducers. An open measurement 
volume refers to a configuration that excludes any obstructions on the upwind side of the 
transducers of an anemometer. This aspect of anemometer design is evident in the 
photograph in Figure (3-1). The open measurement volume concept was not so important 
in the BASS design because the mean flow in the ocean rarely exceeds one meter per 
second. Additionally, the BASS vorticity meter is often used in oscillating flows where it 
is not possible to define a preferred direction. Therefore, the BASS used a center post to 
support its vortical paths. This aspect of the BASS design is evident in Figure (3-2). 

Due to the large mean flows that occur in the atmosphere, it may be necessary to 


create an open measurement volume for the final design of the atmospheric vorticity 
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Figure 3-1. Picture of open measurement volume of an anemometer 
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Figure 3-2. BASS instrument and close view of four transducers per pod 
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meter. The open measurement volume design should provide air that is minimally affected 
by instrument generated wakes for angles of attack that are less than +90°. A possible 


configuration that would permit the open measurement volume to exist in the atmospheric 


Side View 


13.0 cm 15.0 cm 
a 


The Top Struts are thicker in the drawing only. 
This makes the top look different from the bottom. 


Top View 
15.0 cm \ 75 
15.0 cm 15.0 cm 
= 


Figure 3-3. A possible atmospheric design 


vorticity meter is presented in Figure (3-3). Figure (3-3) represents the basic requirements 
for atmospheric use, but the issues of weight and vibration will probably force changes to 
be made. As noted in the top view from Figure (3-3) all acoustic path lengths are 15 cm 


in length. 


There are many advantages to longer path lengths through air from an accuracy point 


of view. However, there are advantages to shorter path lengths from analysis and 
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feasibility points of view. From a feasibility point of view, the received signal strength 
becomes closer to the noise level as path lengths and frequency increase. Once a plane- 


wave sound pulse is in the air, its amplitude is governed by the following equation: 


A =A exp *? (3-1) 


where A is the amplitude of the sound wave, « is the attenuation coefficient in units of per 
centimeter, and x is the distance in centimeters. In fluids, the attenuation coefficient 
increases as the square of the frequency of a sound wave. For dry, CO,-free air at 20°C, 


a has a value determined by the following equation: 
Сю (3-2) 


where f is the frequency of the sound wave in Hertz (Beyer and Letcher, 1969). A value 


for œ in moist air at 30 C is given by: 


&-1.68-107 f7 (3-3) 


(Herzfeld, 1959). Equations (3-2) and (3-3) are based on the frequency squared 
relationship of Stokes-Kirchhoff classical absorption theory. There are factors that cause 
deviations from these curves-of-best-fit to the frequency squared relationship. However, 
the deviant response of absorption due to relative humidity and temperature occurs at 
frequencies that are below our frequencies of interest (Kneser, 1935). In fact, almost all 
of the difference between equations (3-2) and (3-3) at frequencies near 200 kHz could be 
accounted for by the deviant CO, absorption in air as measured by Grossman (Wood, 
1948). The absorption of sound by CO, is maximum at 100 KHz, and it is an order of 
magnitude smaller than the classical absorption at 200 kHz. Although there are 


measurements of absorption coefficients for frequencies near 200 kHz and 1.75 MHz that 
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are higher than equation (3-3), these small factors are unimportant to the analysis that 
follows. Using the more conservative estimate of attenuation given by (3-3), the necessity 
of short path lengths for high frequency sounds in air can be established. Figures (3-4) 
and (3-5) show the attenuation of the amplitude of a plane-wave using equation (3-3) for 


frequencies of 200 kHz and 1.75 MHz. 
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Figure 3-4. 200 kHz attenuation Figure 3-5. 1.75 MHz attenuation 


A straightforward way to examine the attenuation of the signal is to discuss it in terms 
of the extinction distance. The extinction distance is given by the inverse of the 
attenuation coefficient, and it represents the distance at which the amplitude has been 
reduced by 1/е, or a 63% reduction, by absorption. The BASS instrumentation requires 
high frequencies to take advantage of its accuracy. The BASS transducers use 1.75 MHz 
as their transmit frequency. The extinction distance of this frequency in fresh water can be 
more than 3800 cm. (Bergmann, 1938). Using (3-3), the extinction distance for this 
frequency in air is 1.94 cm, as shown in Figure (3-5). Over a 15.0 cm path through moist 


air, a 1.75 MHz sound wave could be expected to lose 99.96% of its amplitude. It is 
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important to note that this expected loss doesn’t include the coupling problem of a 
transducer to air or the spreading losses. 

Acoustic anemometers commonly use frequencies near 200 kHz. Using (3-3), the 
extinction distance for 200 kHz in air is 148.8 cm. Over a 15.0 cm path through moist air, 
a 200 kHz sound wave could be expected to lose 9.59% of its amplitude, as shown in 
Figure (3-4). Clearly, the 200 kHz attenuation is acceptable, but the 1.75 MHz 
attenuation is not acceptable over a 15.0 cm path length. Therefore, it is necessary to 
compromise between acoustic frequency and acoustic path length in air. Assuming the 
signal frequency maintains a value near 200 kHz, the above mentioned 15.0 cm path 
length could be expanded to even greater lengths. 

From an analysis point of view there can be increased difficulties generated by longer 
path lengths. One of the more obvious concerns is vibration of the struts and the frame. 
As the path lengths increase, the strut lengths and frame height would require more length 
if the model shown in Figure (3-3) were used. As the lengths increase, the effective 
stiffness of the members would probably be decreased. This decrease in stiffness would 
eventually lead to vibration that would alter the path lengths. Because the vibrational 
modes would be generated by random forcing functions, it would be extremely difficult to 
remove all of the vibrational-mode errors in the raw data. 

Secondly, if the points raised in the introduction are revisited, another difficulty with 
longer path lengths comes to light. A major objective of this vorticity meter is to analyze 
boundary layers close to the sea surface. For this purpose, a small measurement distance 
between horizontal paths is optimal. If the path lengths were made large enough in the 
vorticity meter, considerably less expensive cup anemometers would provide a more 
accurate measurement of the shear in the boundary layer. The purpose of going to 
acoustic measurements would be lost. 

Assuming that the measurement volume is open (as described above), the percentage 
of an acoustic path that would be affected by transducer wakes would be reduced by 
longer path lengths. In an open measurement volume, the acceleration of flow around the 
transducers is interpreted as a faster or slower velocity in the two portions of the acoustic 


paths that are affected. These local increases in velocity cause a small change in the 
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averaged acoustic velocity over the path. The simplest and best way to correct for this 
acceleration problem is to perform wind tunnel tests and correct the results empirically 
(e.g., Grelle and Lindroth, 1994). Of course, path lengths that are long enough would 
make the empirical corrections negligible. In addition, the small errors in the timing 
circuitry would represent a smaller portion of the total time that an acoustic signal would 
require to transit the acoustic path length. Despite these advantages, the acoustic path 
length requirements and limitations prevent a lengthening of acoustic path lengths in the 
initial instrument design. Therefore, a signal of approximately 200 kHz over a 15.0 cm 
path length is an adequate compromise. However, it is unlikely that the acoustic vorticity 
meter will be used in only one acoustic path length as it evolves. 

The empirical corrections to velocity measurements from a wind tunnel test are most 
useful if the instrument is as symmetric as possible. If the wake generation effects on 
velocity paths were not symmetric, the instrument would require calibrations at various 
angles to the mean wind. This would require that an estimated mean wind be backed out 
of the measurements before corrections could be applied. This is the reason that the six 
transducer pods in Figure (3-3) have a spherical shape. Each of the six transducer pods 
shown in Figure (3-3) would contain four transducers in a configuration that would be 


similar to the pod configuration in BASS. 


TRANSDUCER DESIGN 

The first difficulty that occurs when attempting to insonify air is an impedance 
mismatch between the transducer and the air. This difficulty occurs because the density of 
piezoelectric ceramics and crystals is much higher than that of air. The effective 
impedance mismatch can be substantially reduced by the use of piezopolymers or the use 
of multiple matching layers covering a ceramic or crystal (Wallace, 1992). Other methods 
use piezoceramic embedded in epoxy to decrease the impedance mismatch between air and 
a transducer (Hayward et al., 1992). 

Our group attempted to use the experience of atmospheric transducer design by 
searching for off-the-shelf transducers that would work with the BASS electronics. The 


first two transducers tested were simple piezoelectric discs that resonated at 1.75 MHZ 
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and 875 kHz. These discs were tested because they were already available at WHOI, such 
that these discs would have represented substantial savings if they showed any promise. 
These discs were tested using a storage oscilloscope with 10x probes. Not surprisingly, 
given the above discussion, the results for these discs were negative. However, it should 
be noted that no amplification or layer matching was used. 

The next attempt used piezoceramic cylinders that were designed for a resonance of 
approximately 200 kHz. The resonance frequencies of the ten cylinders that were ordered 
were found to vary between 240 kHz and 310 kHz. In an attempt to decrease the 
impedance mismatch between the cylinders and air, RTV was used as a filling material for 
the cylinders. Once the cylinders were filled with RTV, the resonance frequencies were 
found to vary from 250 kHz to 330 kHz. Any two cylinders were separated by a distance 
of 0.15 m, and one transmitted under a square wave voltage with an amplitude of 20 volts 
while the other received. Once an amplifier with a gain of 480 was developed, a signal of 
approximately 30 mV was visible on the storage oscilloscope. However, these 
transducers were found to be unacceptable because of their large variance in resonance 
frequency and because of their poor coupling with the air. 

The final transducer design tested was the E-188/220 by Massa Products. It has an 
estimated 220 kHz resonance frequency, and it takes advantage of impedance matching 
materials. The E-188 consists of a piezoceramic disc with a quarter wavelength silicon 
matching layer. This transducer design was much more dependable in its proximity to the 
expected resonance frequency, and the E-188 was found to work well in air. Without 
amplification, an E-188 transmitting under a square wave of frequency 220 kHz and 
amplitude one volt can generate up to a 1.5 mV signal in another E-188 that is 0.15 
meters away. These results clearly indicate that the design of the E-188 allows a coupling 
to air that is orders of magnitude more efficient than the coupling of the piezoceramic 


cylinders. 


RINGING ANALYSIS 
Because the estimated length of separation of the transducers was expected to be 0.15 


m in the test design, the continued oscillation, or ringing, of the transducers after the 15- 
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cycle sine wave pulse was a concern. In order to test the ringing of the E-188, the 
transducers were modeled as a transfer function and were separately tested against an 
oscilloscope. The initial modeling of the transfer function was performed on the E-188's 
with a Hewlett Packard 4195A Network Analyzer. In our tests, the E-188's were 
subjected to a multi-cycle sine wave, and the results were recorded on a storage 
oscilloscope. 

The Hewlett Packard Network Analyzer 1s capable of modeling crystal oscillators using 


the following equivalent circuit: 








Figure 3-6. HP ringing model 


where L is an inductor, R is a resistor, and C is a capacitor with subscripts denoting 
different values. This equivalent circuit is capable of generating the following Laplace 


transfer function: 
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(3-4) 


However, the addition of a quarter wavelength silicon layer to the piezoceramic disk 
creates a situation where two separate systems interact. This overlapping of systems can 


be modeled by the following parallel equivalent circuit: 


23 











Figure 3-7. Ringing model 2 


where C, ‘s are capacitors, L,’s are inductors, and R,’s are resistors. The overlapping of 


two crystal oscillators generates the following Laplace transfer function: 
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This result can be rewritten by creating lıke denominators and removing them from the 


function. This yields the following Laplace transfer function: 
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This form of the transfer function allows the power of the numerator and denominator to 
be easily seen. Since the power of the denominator is larger than the power of the 


numerator, a first approximation suggests that the system should be stable. 


The Hewlett Packard Network Analyzer was able to produce the following values with 


its equivalent circuit analysis option for crystal oscillators like Figure (3-6): 


R= 45.5342 Q 

L = 0.980964 mH 
C,= 597.691 pF 
Cos ВА ПЕ 


Using these results in equation (3-4) produces the following bode plot: 


Hewlett Packard Model 


Mag. Impedance 





Frequency (rad/sec) 


Phase deg 





Frequency (rad/sec) 


Figure 3-8. Bode plot of Hewlett Packard model 
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The ‘Hewlett Packard Model’ plot has some similarity to the manufacture’s specifications 
and the measurements displayed by the network analyzer, but it is missing a second peak 
near 220 kHz. By modifying the values provided by the network analyzer, it is possible to 
produce a bode plot that resembles the manufacturer’s specifications and the actual 
measurements displayed by the network analyzer. The modified values for use in equation 
(3-6) are as follows: 

R= 45.5342 Q 

L = 0.980964 mH 

C,= 597.691 pF 

C,= 1.23217 nF 


К„= 59 С 
L,= 0.84 mH 
= 975 pF 
a 2 nF 


The first four values are ıdentical to the four measured values generated by the network 
analyzer for a crystal oscillator. However, the use of these eight values in the transform 


function of the overlapping system produces the bode plot shown in the following figures: 


Transducer Model 
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Frequency (rad/sec) 


Figure 3-9. Bode plot for transducer model 2 
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This “Transducer Model’ bode plot closely resembles the manufacture’s specifications for 
the E-188/220. By placing the equivalent circuit in Figure (3-7) under a 14.2449-cycle, 
220 kHz sine wave of amplitude one volt, the following plots of the transducer input and 


transducer response were generated: 


Transducer Input 





о 0.5 1 1.5 2 25 3 
Time (sec) x 107% 


Figure 3-10. Assumed input to an E-188/220 transducer 


Transducer Ringing 
150 


Response to Signal 





Time (sec) х 10-* 


Figure 3-11. Hypothetical output of transducer using assumed input 
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These analytıcal results show a transducer response that stabilizes at some voltage offset. 
The voltage offset changes as the number of frequency cycles changes. For example, the 
use of a 14-cycle, 220 kHz sine wave would show a more negative voltage offset in Figure 
(3-11). The use of a 14.5-cycle, 220 kHz sine wave would show a small positive offset 
that would be almost equal in magnitude to the offset shown in Figure (3-11). 

These results do not make physical sense. However, the ringing response about the 
stabilization voltage looks the same in all of the cycle lengths used in our tests. Therefore, 
there may be some value in the analysis. The analytical results are dubious by any 
standard because the conditioning number for the values used in equation (3-6) to 
generate the transducer response was almost 4X10. Therefore, it was necessary to look 
at the physical characteristics of the transducer ringing after a multi-cycle input. More 
information on the analytical portion of this analysis is presented in Appendix A. 

Because the ringing of the transducers was still a concern for measurements using 
modified BASS instrumentation, a physical ringing test was set up in such a manner that it 
would determine whether ringing was a concern or not aconcern. Two transducers were 
set up at a separation distance of approximately 15 cm. The transmitting transducer was 
connected to a circuit with the effect of a one shot, i.e., the length of time between 
transmissions was substantially longer than the length of transmission. This ‘one shot’ was 
tuned for a multi-cycle sine wave of approximately 220 kHz, and a schematic of its design 
is shown in Figure (3-12). 

The output of the circuit was tuned for a multi-cycle sine wave at 220 kHz every 
1/1000" of a second. The length of the multi-cycle sine wave was almost ten times shorter 
than the length of time between pulses, which provided an effective one shot to the 
signaling transducer. The output of the receiving transducer was measured on a LeCroy 
Digital Oscilloscope. The Lecroy Digital Oscilloscope was connected to 10X probes that 
provided the only amplification of the transmitted and received signals. These signals are 
displayed in Figure (3-13) and Figure (3-14), respectively. The transmitting signal had a 


duration of 75 ps, and although there is some physical ringing, it rapidly falls near the 
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noise level. In air, a sound wave will take about 440 ts to traverse a 15.0 cm pathlength. 


Therefore, the ringing of the transducers did not appear to be a concern. 
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Figure 3-12. Circuit used to generate a short burst of 220 kHz waves 
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Figure 3-13. Transmitted Sıgnal through an E-188/220 
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Figure 3-14. Received Signal using an E-188/220 
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TIMING AND BURST GENERATION 

Due to the different group velocities of sound waves in air and water, the timing for a 
pulse in air over a 15.0 cm path length must be different from the timing of a pulse in 
water over a 15.0 cm path length. In addition to the timing changes, the resonant 
frequency had to be changed from 1.75 MHz to 220 kHz. In the BASS instrumentation, 
the timing and resonant frequency generation is performed within the circuitry of the 
timing and burst generation board, and its original form is displayed in Appendix B. 
Although the alterations to the board required several iterations, only the final changes are 
detailed in this section. 

In order to change the resonance of the board, a resonant LC circuit had to be changed 
to meet the requirements of 220 kHz transducers. This LC circuit is represented in 
Appendix B by capacitor C, and inductor L,. Resonance in the LC circuit follows the 


following equation: 


1 


€ | 
2x LC, iu 


where / is the frequency in Hertz. The divide by two behavior of the circuit creates an 
additional factor of approximately 2.26. Therefore, the new LC circuit should combine an 
L, and C, that are resonant near 495 kHz. Because the availability of capacitors exceeds 
that of inductors, C, was replaced with a 3600 pF capacitor. This produces a resonant 
frequency of slightly more than 218 kHz going to the E-188/220 transducers. 

The timing of the circuitry also had to be changed. During these modifications, we 
were able to take advantage of a previously unused clock that was available on the board 
(represented by U9B on the schematic in Appendix B). In the air-side circuitry, that clock 
was used to double the length of each clock step: 1.e., one clock cycle was changed from 
6.5 ps to 13 us in length. This change allowed the new timing board to compensate for 
the generally higher velocities found in wind measurements as compared to current 
measurements in water. The counters U, and U, on the timing and burst generation board 


were used to alter the settings for oscillator control, start transmit time, transmit disable, 


start A/D, and reset. 
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The oscillator control is turned on first so that any capacitance is charged in the 
resonant part of the circuit before a signal is sent to a pair of transducers. Next, the start 
transmit goes high to set the circuitry such that it will transmit the resonant pulse when the 
clock resets are turned off. When transmit disable goes low immediately after start 
transmit goes high, the clock resets are turned off, and pulse transmission begins. This 
pulse is generated by the resonant part of the circuit, and it is transmitted to a pair of 
transducers as long as start transmit is high. Start transmit eventually goes low, and the 
transducers are no longer excited by the circuitry. After a short time, the unclamp 
receiver goes low, which allows the same transducers to act as receivers and to 
communicate their received signals to a set of counters. When fourteen zero-crossings 
have occurred at one counter due to the cycling of its transducer at the resonant 
frequency, the counter switches a constant current source to an integrating capacitor. 
Assuming there is some flow along the acoustic path, the other counter will count fourteen 
zero-crossings from its transducer at a slightly later time, and that counter will switch a 
constant current source to another integrating capacitor. After some voltage has built on 
the two integrating capacitors, the difference between the values is recorded on the A/D 
board when start A/D goes high. The system is then reset, and another path is measured. 


The cycling of these settings is represented in the following figure: 
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Figure 3-15. Cycling of board 
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The cycles in Figure (3-15) were created by using the following U, and U, connections: 
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Figure 3-16. Alterations to timing and burst generation counters 


The BASS instrumentation contains a Tattletale which is used to make an equivalent 
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time difference based on the voltage difference that is recorded on the A/D board. The 
Tattletale is programmed through a language called TTBasic, and a program similar to the 
ones used by BASS instrumentation is provided in Appendix D. Once an equivalent time 
difference is calculated, it is used to calculate a velocity difference. The measurement of 
the velocity difference can be improved by iterative processes. If these iterations were 
done by the TTBasic program, the Tattletale would be in non-sleep modes for longer 
periods, and this could cause power requirement problems for longer, battery powered 
deployments. If a long deployment only returned voltage difference or time difference 
data, the processing of twelve paths sampled at 10 Hz for 100 days could require more 
processing time that a researcher would like to spend. The error analysis in Chapter 4 
assumes the use of the simple equation for velocity that does not require iteration and is 


already used by BASS. 
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CHAPTER 4 


ERROR ANALYSIS 

Typical commercially available sonic anemometers compute the time it takes for a 
sonic pulse to travel between two transducers. The time of flight for the pulse is altered if 
the wind has a velocity component along the path between the transducers. For two 


pulses traveling in opposite directions, the times of flight are given by: 





/ EE 4 
mem (4-1а) 
E d 

N (4-1b) 


where t, and t, are travel times in opposite directions, d is the acoustic path length, c is the 
local speed of sound in the medium, and v is the spatially averaged velocity component 
along the acoustic path. The wind velocity can then be determined from the difference of 


the reciprocal times by : 


т (4-2) 


The advantage of this method is the removal of the speed of sound from the equation. 
BASS does not record the travel times separately, it uses the difference between travel 
times to determine the flow component along an acoustic path. The time difference can be 
used in the following equation to calculate velocity: 
d (al 


=- —+ + 4-3 
на те (4-3) 
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where the correct root is easy to determine since one gives a velocity greater than the 
speed of sound. However, the accuracy of this equation obviously rest with an accurate 
measurement of the speed of sound. Equations (4-1a) and (4-1b) can be combined to 


yield an expression for the difference in travel times as: 
2 24 
К ы = А Sc) 
2 


Some BASS TTBasic programs use only the first term in equation (4-4), such that the 


velocity is derived from: 


„ЄЛ 
у= 
24 





(4-5) 


The error caused by the exclusion of the higher order terms for measurements performed 
in water 1s approximately 1 ppm (Williams et al., 1987). Because the ratio of air velocity 
to speed of sound 1s generally larger in air, we expect this error to be substantially larger 
as well. The purpose of the following sections is to demonstrate where the more 


important errors will or can exist so that they may be eliminated if necessary. 


ERRORS FROM HIGHER ORDER TERMS 

The errors from higher order terms in equation (4-4) are more substantial in air than in 
water for two reasons. First, the average velocities seen in air are higher than those seen 
in the water. Second, the speed of sound in air is slower than the speed of sound in water 
by a factor of almost five. Therefore, the errors generated by the exclusion of the higher 
order terms in equation (4-4) require evaluation in order to determine whether (4-5) is a 
valid approximation in the atmosphere. The results of a comparison made at a virtual 


temperature of 273 K with an acoustic path length of 0.15 m are presented below: 
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| Velocity | Time Diff, Time Diff, | % Error | 
























| 3Terms| 1Term]| A 
(m/sec) | (msec) (msec) | ЖҰ 
. 1.00000 0.00273 0.00273 | 0.00091 
0.0136: 9082202 

ОВЕ _| 0.0908 

_ 0.04090 .| 0.20451 

05474 0.05454 | 03635 

0.06856 0.06817 _ ^| 05680 
0.08248 0.08180 .— — | 0.8179 


Table 4-1. Error ın higher order terms 


It is an obvious result from the table that errors due to the omission of the higher order 
terms should not exceed one percent of the true value of the time difference. Since the 
time difference and velocity are linearly related, the same can be said for our velocity 
estimates. Thus, the higher order terms in equation (4-4) will be neglected during this first 
stage in the development of the atmospheric vorticity meter. If more accuracy than that 
provided by equation (4-5) is required, a combination of temperature and/or hygrometer 
measurements could be used to deliver an accurate speed of sound for use in equation (4- 
3). Equation (4-5) is used in the error analysis because it is used by BASS, and it provides 


more computational efficiency than equations (4-3) or (4-4). 


ERRORS DUE TO VIRTUAL TEMPERATURE CHANGES 
The speed of sound in air is a function of both temperature and humidity. These effects 
are most easily related to the speed of sound using the sonic virtual temperature in the 


following equations: 


c= (4037, (4-6) 
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Г, =1(1+0.514) (4-7) 


where c is the group speed of sound waves, 7; is the sonic virtual temperature, T is 
temperature and q is specific humidity (Schotanus, 1983). This quantity differs slightly 


from the classically derived virtual temperature given by: 


f =/(1+0.61q) (4-8) 


where 7, is the virtual temperature. Virtual temperature is the temperature that dry air 
would need to have a density equal to moist air at the same pressure (Stull, 1988). 
Because water vapor is less dense than diatomic oxygen and nitrogen, a packet of moist 
air would be less dense than a packet of dry air at the same temperature and pressure. 

Although the ability to neglected the higher order terms in equation (4-4) has been 
shown, the calculations that were involved assumed perfect measurements of virtual 
temperature. Accurate measurement of humidity is often difficult in the marine 
environment. Therefore, it is of interest to see what the exclusion of this correction would 
do to our measurements. Assuming that a proper measurement of temperature is provided 
at 300 K, the following tables show the errors that might occur in equation (4-5) for 


relative humidities between 30% and 100%: 










Relative | Velocity | Velocity | % Error 















Humidity | Temp Temp & 
(%) Humidity s 
30.00 5.00 5.02 





40.00|. 500 50) | 
__ 50.00 5.00 5.03 |_ | 0.55 
60.00 5.00 BUE 0.67 


70.00] | 500| 504 0.78 
80.00 5.00 5.04 


90.00 5200 5.05 1.00 


100.001 — 5.00 5.06 ER 


Table 4-2. Humidity error when measuring 5 m/s 
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| Relative Velocity | Velocity | E Error 
IAS | | 
EU PM | Humidity | | | 
_ 30.00] 1000] 1003 ES 
40.00; 1000 1004 | 045. 
_ 50.00 1000 
| . 60.00 10.00 








___ "000 
80.00 


_ 90.00 
| 100.00 


10.00 
10.00 
10.00 
_ 10.00 























Table 4-3. Humidity error when measuring 10 m/s 


Once again, the errors produced are much less than 2%, even for extreme cases. 
Therefore, the use of a separate hygrometer is not necessary for reasonable accuracy in the 
acoustic velocity measurements. In fact, it is worth noting that the errors resulting from 
the exclusion of higher order terms in (4-5) and specific humidity in (4-6) tend to offset 


each other. 


ERRORS DUE TO FLUCTUATING TEMPERATURE 

Since the need for a separate hygrometer and the use of higher order terms in equation 
(4-4) have been eliminated, the next step 1s a test for the need of an atmospheric 
temperature sensor. The vorticity meter should be capable of operating in a broad range 
of temperatures. Therefore, an operating range of 273- 310 K will be assumed. Ifa user 
were trying to pick a temperature that would minimize errors throughout this range, an 
fixed temperature of 293 K might be used. The following table compares the errors 
generated by this temperature assumption when a 0.15 m acoustic path is subjected to the 
extreme temperatures of its range. 

In Table (4-4), the time differences are first computed by inverting equation (4-5) (as 
if the modified BASS instrument was operating in air) and using the velocity that is shown 
in the table to calculate the speed of sound at a temperature of 273 K through equation (4- 
6). The following column shows the velocity calculated using the same time difference 


with an assumed temperature of 292 K. The percent error is presented in the next column. 
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. Velocity. Measured Velocity | % Error | ! Measured Velocity | % Error 
(m/s) | Time Diff, w/ 292K. in 292 K | | | Time M w/ 292 K { in 292 K 
в... 273К (ms / Assumed e 2/73. 310К (ms, Assumed wi T= 310 
1.00000, 0.00273| 1.06960 6.959711, _ 1740; 00240 0.94193 5.80658 
5 00000. 0.01363! 5.34798| 6.95971 | E 4.70968 5.80634. 
10.00000 0.02727| 10 p 6. seri 0.02401| 9.41936 
15.00000 0.04090! 16.04396 6.95971 0.03602} 14.12903 
20. 60000. 0. 05454 21.39194 6.95971 | |. | 0.04803] 18.83871 
7| 26.73993| 6.95971 _ | 0.060031 23.54838 5.80647 
0.08180 32.08791| 6.95971 0.072041 28.25807 _ 5.80644 





Table 4-4 Assumed temperature error 


The next three columns repeat the process for a measurement taken at a temperature of 
310K. For the acoustic paths to meet the accuracy requirement, it is evident that a 
temperature sensor of some form will be necessary. This temperature sensor would be 


used to adjust the speed of sound estimate generated by equation (4-6). 


ERRORS IN THE MEASURED TIME DIFFERENCE 

Throughout the error analysis, tt was assumed that the time difference as measured by 
the BASS instrumentation was without substantial error. For BASS time difference 
measurements, the accuracy is approximately 40 ps (Trivett, 1991). The extreme accuracy 
of this measurement is generated by the cancellation of noise generated errors within the 
BASS instrument when the voltages from the integrating capacitors are subtracted. The 
40 ps error in time differences is several orders of magnitude smaller than the time 
difference generated by a 1 m/s mean wind. Therefore, this error is not considered. 

BASS instrumentation has shown an accuracy of 7 ns for time of flight measurements 
(Trivett, 1991). As shown by equation (4-2), the use of times of flight instead of time 
differences would alleviate the need for additional temperature and/or humidity sensors 
when taking velocity measurements. A 7 ns time error would equate to a velocity error of 
less than 0.6 cm/sec, i.e., it would provide a much better estimate than the errors shown in 
the rest of this chapter. However, the time of flight technology is being put into BASS 


instruments at the time that this thesis is being written, and it is not immediately available. 
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СНАРТЕК 5 


WINDTUNNEL 

A windtunnel was used to calibrate and test the velocity measurements of one axis of 
the atmospheric vorticity meter. In order to test one axis of the vorticity meter properly, it 
was necessary to understand the characteristics of the windtunnel that was used. The 
velocities at three cross-sections, with twenty-eight experimental samples per cross- 
section per wind speed, were measured with a new pitot tube. The cross-sections were 
taken at distances of 2.5 cm, 71.8 cm, and 98.4 cm behind the front, leftmost, window of 
the windtunnel shown in the photographs. The photographs in Figure (5-1) display the 
experimental section and shape of the windtunnel. The fan is used to puil air, rather than 
push air, through the tunnel in order to generate a more laminar flow field inside the 
experimental section. The experimental section has a height of 61.0 cm and a width of 
91.4cm. The air intake at the front of the tunnel has a height of 1.05 m and a width of 
1.89 m. 

The wind-speeds at the center of the second cross-section were used as the normalizing 
velocity for each of the three velocities used in the calibration. The three, normalized 
velocities were 2.35 m/s, 4.37 m/s, and 8.81 m/s. Using the pitot tube, an experimental 
sample was taken approximately every 10.2 cm in the horizontal and every 7.6 cm in the 
vertical for each cross-section. The experimental samples avoided the top and bottom 7.6 
cm as well as the first 10.2 cm on each side of the windtunnel. 

The windtunnel had no record of calibrations except the one that occurred when the 
windtunnel was first built. That calibration was thirty years old. Therefore, there was 
some concern that the characteristics of the windtunnel had changed. When the new 
calibration experiments were performed, the windtunnel displayed the same characteristics 
that it had displayed thirty years ago. Based on this thirty-year stability, 1t was assumed 
that the windtunnel results were stable for the few weeks of its use during this experiment. 
A contour plot of the windtunnel’s second experimental cross-section at a normalized 


velocity of 2.35 m/s is provided in the Figure (5-2). 


41 








Figure 5-1. Photographs of the windtunnel 
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Figure 5-2. Contour plot of second cross-section 


Appendix C contains the Basic and C programming files that were written for use with 
the windtunnel. The parts of the programs which refer to CIO8 are referring to files for 
the A/D board. The A/D board files came with the A/D board, and they are not included 
for this reason. In addition, a portion of the CPU’s RAM was defined as the D drive when 
Windram.bas was running. This was done in the config.sys file using ramdrive.sys. The 
use of RAM as a dummy drive allowed more rapid input and output of data to and from 


the file Diff.tmp. 
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СНАРТЕК 6 


RESULTS AND CONCLUSIONS 

This thesis was a test of the feasibility of converting BASS instrumentation for the 
ocean into atmospheric instrumentation. There were several changes to the circuitry that 
were expected for the transfer of technology to be successful. First, the resonant 
frequency of the transducers needed to be reduced because the absorption of high 
frequency sound in air is greater than that for water. Second, the timing of the circuitry 
needed to be adjusted in order to account for the slower group speed of sound waves in 
air in comparison to water. Third, the resonant frequency of the LC circuit used to excite 
the transducers had to be lowered to the resonant frequency of the new transducers. 
Using these few changes, it was assumed that a working instrument for atmospheric use 
could be developed. 

After an extensive search for a set of transducers that were resonant at a frequency 
near 200 kHz, the E-188/220's by Massa Products were selected. These transducers 
provided adequate coupling to air due to their 1/4 wavelength silicon matching layer. The 
optimal signal frequency of these transducers occurred near 220 kHz, and the four 
transducers of the initial order compared nicely in their response to an HP4195A 
Spectrum Analyzer. The E-188/220's were then subjected to a ringing analysis. The 
results shown in Figures (3-13) and (3-14) show that the ringing of the transducers for 
path lengths of 15.0 cm would not be a concern. 

An error analysis concerning the use of the current BASS equation for velocity and the 
accuracy of the time difference measurement was also conducted. This analysis showed 
that the equation used for analyzing BASS velocity data would be of sufficient accuracy 
for some tests. However, the use of more terms in equation (4-4) or equation (4-3) in the 
TTBasic program could prove to be a quick way to increase accuracy for air 
measurements. The error in time difference measurements for BASS instrumentation is so 
small that it is not worth consideration at this point. 

Finally, a windtunnel was prepared to test the accuracy of the velocity measurements 


against the measurements made by a calibrated pitot tube. A QBasic code which recorded 
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the time of a measurement and the temperature in °C, barometric pressure in mbars, and 

velocity in m/s of that measurement was developed. Additionally, the windtunnel showed 
very long term stability in its calibration measurements. The measurements of velocity (or 
simply a voltage difference) at a particular time for BASS instrumentation are recorded on 


its internal Tattletale using TTBasic. 


INITIAL BENCH TESTS 

When the modified timing and burst generation board was inserted into the remaining 
BASS circuitry, there were two problems. First, the almost 60 dB losses of the E- 
188/220's over a 15.0 centimeter path length were higher than those of the BASS 
transducers through water. In this case, the loss was mainly due to the impedance 
mismatch and spreading loss that resulted in a low received signal (rather than plane-wave 
attenuation of the signal). Secondly, there was noise in the receiver that interfered with 
the received signal. This noise was occurring in the same frequency band as the received 
signal. These problems had to be fixed before a working device could be constructed. 

The problem that could be solved was that of signal loss. The signal in the BASS 
instrumentation limited the input to a transducer to four volts. The signal was increased to 
twelve volts and a step-up transformer was used to achieve a twenty-nine volt transmit 
signal. The use of a twenty-nine volt signal allowed the receiving transducer to receive a 
signal near 300 mV. This signal was large enough to be detected in the zero-crossing 
counting circuitry. Therefore, simply generating a larger signal solved the first problem. 

The second problem was more complex. The circuitry ringing was unexpected because 
the ringing analysis had shown rapid reduction in signal strength after the transducers were 
removed from the signal. Therefore, we suspected that this problem was due to ringing 
inside the BASS instrumentation. In order to be sure, we double checked the ringing 
analysis for the transducers and obtained results that were very similar to the previous set 
of data. Therefore, there was some other problem. 

Because the LeCroy analyzer was having no problem detecting a signal, the BASS 
circuitry had to have a different nature that allowed the ringing to occur. The first 


attempts to solve the problem involved matching the impedance for both the transmit and 
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receive sides of the circuit at a lower impedance. The impedances needed to be matched 
to avoid phase shifts, i.e., when voltages are zero, the phase shifts between voltages and 
currents are zero (Brown and Lawson). The impedance was reduced in order to get a 
larger signal into the air. This produced slight improvements, but they were not enough to 
detect a received signal. Another attempt involved adding circuitry between the 
transducers and the existing BASS circuitry that was designed to damp the receiver noise. 


The circuit and its relationship to a transducer are shown in Figure (6-1). This circuit 
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Figure 6-1. Damping circuit 
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worked nicely on a temporary (bread) board, but when it was incorporated in the BASS 
instrumentation, a received signal could not be detected above the circuit noise. Due to 
the time constraints involved, it was proving to be unlikely that an actual measurement 
would be made. 

The next possible solution is to increase the impedance going into the transducer in 
order to create more damping in the circuitry. Of course, this introduces the possibility of 
phase shifts between voltage and current which would introduce errors into measurements 
of the time difference. These errors could be calibrated out of the measurement. Since 
the transducers were not tested for fluctuations due to temperature shifts, the calibration 


might need to be conditioned on temperature (Brown and Lawson). 


FUTURE TUNNEL AND FIELD TESTS 

The first test that will need to be performed with BASS instrumentation is a one path 
velocity measurement in the windtunnel. When the noise in the circuitry 1s small enough 
that a received signal can be detected, the rest of the equipment for these tests will be 
waiting. The first result to look for will be voltage outputs (as recorded in the Tattletale) 
that are proportional to the different wind speeds. These windtunnel tests should not 
exceed 2.4 m/s until corrections are made to linearize the constant current source response 
with faster transistors (Williams, 1995). If the changes in velocity prove to be 
proportional to the changes in voltage recorded by BASS, the next step will be to 
determine the cosine response of the acoustic path axis to the direction of velocity. When 
the cosine response measurements are made, the effects of flow distortion by the 
transducers should be apparent. 

Errors in the velocity measurements due to flow distortion should show proportionality 


to the following equations (Thwaites, 1995): 


Ud | 
ES os0 Potential Flow (6-1) 
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3-D Sphere Wake (6-2) 


where d is the diameter of the object in the flow, / is the acoustic path length, and L is the 
distance downstream of the object. However, it has been noted that it is better to calibrate 
structures in a windtunnel than to attempt calibration using theory (Wyngaard, 1981). If 
left uncalibrated, these flow distortion errors would make accurate shear measurements 
impossible. Even calibrated acoustic instruments can have errors as large as 45 cm/sec for 
ten second averages of wind speed (Solent Research, 4.1). However, these errors are 
generated by instruments that do not have the accuracy of time measurements that BASS 
does. Additionally, Thwaites (1995) has shown that wake generated errors tend to cancel 
using a four-path design. Ultimately, the accuracy of a four-path axis calibration will 
determine if a velocity shear measurement can be made with 15.0 cm acoustic paths. 

If 15.0 cm path lengths are insufficient to provide shear measurements, the transducers 
could be separated by longer acoustic paths. Equations (6-1) and (6-2) show that the 
acoustic path length is inversely proportional to the expected flow distortion errors. In 
addition, the differences in velocities created by shear would be greater over longer paths. 
However, the error created by not knowing the zero crossings precisely would increase 


with path length according to the following (Helstrom, 1975): 


Voltage Noise 


- 2llf «(Signal Voltage) Sa 


due to the attenuation and spreading of the signal over longer path lengths. 

If the 15.0 cm paths аге sufficient to measure shear, a four-path axis would then be 
deployed for field tests. The horizontal acoustic paths would need to be aligned with the 
mean wind, and they would be used to measure shear. The measured difference between 
horizontal velocity in the top and bottom paths could be used in equation (2-3) to 


determine values of u.. These values could be compared to the u. values generated by a 
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sonic anemometer deployed nearby. 

Eventually, a three-axis vorticity meter could be deployed. It could provide 
simultaneous direct measurements of shear and u. in the same measurement volume. This 
would allow an improvement in the estimates of von Karman’s constant. The three-axis 
vorticity meter could be deployed under various conditions to determine the most likely 


value for von Karman’s constant. 


49 





References 


L. Bergmann, Ultrasonics: And Their Scientific and Technical Applications, John Wiley 
&Sons, New York, 1938. 

R. T. Beyer and S. V. Letcher, Physical Ultrasonics, Academic Press, New York, 1969. 

N. L. Brown and K. D. Lawson, A High Precision Acoustic Current Sensor, Neil Brown 
Instrument Systems, Cataumet, MA. 

A. Grelle and A. Lindroth, “Flow distortion by a solent sonic anemometer: wind tunnel 
calibration and its assessment for flux measurements over forest and field," J. Atmos. 
Oceanic Techn., 11, 1529-1542, 1994. 

G. Hayward, A. Gachagan, R. Hamilton, D. A. Hutchins, and W. M. D. Wright, 
"Ceramic-epoxy composite transducers for non-contacting ultrasonic applications," in 
New Developments in Ultrasonic Transducers and Transducer Systems, Proceedings of 
the SPIE conference held 21-22 July 1992, Belingham, WA, pp.49-56. 

K. F. Herzfeld, Absorption and Dispersion of Ultrasonic Waves, Academic Press, New 
York, 1959. 

J. R. Holton, An Introduction to Dynamic Meteorolog y, Academic Press: International 
Geophysics Series, New York, 1979. 

H. O. Kneser, “Molekulare schallabsorption in gasen”, Z. Techn. Phys., 16, 213, 1935. 

P. Schotanus, F. T. M. Nieuwstadt, and H. A. R. de Bruin, "Temperature measurement 
with a sonic anemometer and its application to heat and moisture fluxes", Bound.-Layer 
Meteor., 26, 81-93, 1983. 

W. A. Smith, "New applications in ultrasonic transducers emerging from innovations in 
piezoelectric materials,” in New Developments in Ultrasonic Transducers and 
Transducer Systems, Proceedings of the SPIE conference held 21-22 July 1992, 
Belingham, WA, pp.3-26. 

Solent Research, Ultrasonic Anemometer, Product Spec., Needham, MA, Issue 4.1 

R. B. Stull, An Introduction to Boundary Layer Meteorology, Kluwer Academic Pubs., 
Boston, MA, 1988. 

H. Tennekes and J. L. Lumley, A First Course in Turbulence, MIT Press, Cambridge, 


50 





МА, 1980. 

F. T. Thwaites, Development of an Acoustic Vorticity Meter to Measure Shear in Ocean- 
Boundary Layers, P.H.D. Thesis, M.I.T.-W.H.O.I., 1995. 

Е. T. Thwaites, A. J. Williams 3rd, E. A. Terray, and J. H. Trowbridge, “A family of 
acoustic vorticity meters to measure ocean boundary layer shear,” in the Fifth Working 
Conference on Current Measurement, Proceedings of the IEEE conference held 7-9 
February 1995, St. Petersburg, FL, pp. 193-198. 

D. A. Trivett, Diffuse Flow from Hydrothermal Vents, P.H.D. Thesis, M.I.T.-W.H.O.L., 
1991. 

A. J. Williams 3rd, J. S. Tochko, R. L. Koehler, W. D. Grant, T. F. Gross, and C. V. R. 
Dunn, “Measurement of turbulence in the oceanic bottom boundary layer with an 
acoustic current meter array,” J. Atmos.Oceanic Techn., 2, 312-327, 1987. 

A. J. Williams 3rd, “Linearity and noise in differential travel time acoustic velocity 
measurement’, in the Fifth Working Conference on Current Measurement, Proceedings 
of the IEEE conference held 7-9 February 1995, St. Petersburg, FL, pp. 216-219. 

R. W. Wood, Supersonics: The Science of Inaudible Sounds, Brown University: Charles 
K. Colver Lectures (1937), Providence, RI, 1948. 

J.C. Wyngaard, “The effects of probe-induced flow distortion on atmospheric turbulence 


measurements", J. App. Meteor., 20, 784-794, 1981. 


ЭД 





APPENDIX A 


Transducer Response Graphs 
A-1 Model using Thevenin equivalent circuit for receiving transducer 
A-2 Direct measurements of transducer frequency response from HP4195A 


A-3 Response of the transducer to an input signal using the direct measurements in A-2 


The forms of ringing analysis that were not detailed in the text consisted of modeling 
one transducer and two transducer systems. The one transducer model of the frequency 
measurements from the HP4195A is shown in A-2 and A-3. The one transducer model 
produces more believable results than the one in the text, but its conditioning number is 
approximately 4.7 X 10? after digitization. The two transducer model consisted of a 
Thevinin equivalent circuit that was generated by modeling the piezoceramic resistor 
ringing into the silicon layer resistor and convolving the result. The Matlab programs that 
were used to generate the graphs immediately follow the graphs that they generated. The 
last two programs in this section show the Matlab routines used to create the graphs in the 


Ringing Analysis section of Chapter 3. 
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х 10714 Ringing for Receiving Transducer Resistor 


Resistor Response 





Time (sec) des 
Figure A-1. Model using Thevinin equivalent circuit for receiving transducer 


To To To To To Yo To To Fo To To To Fo Fo Fo Fo Fo To To To To Fo Fo Fo Fo Vo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo To Lo To To To Vo 
% File saved as ring_eq.m that was used to generate the receiving transducer ringing plot 
clear 

L= 0.980964e-3; 

L2= 0.00084; 

R= 45.5342; 

R2= 35; 

Сі- 5.97691е-10; 

C2= 1.23217e-9; 


C3= 5.75e-10; 

C4= 2e-9; 

% Thevin equivalent about resistor for Silicon layer only 
Jonum= [1 0); 9o ] over 
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Yden= [R2*L2 1 R2/C3]; % 5+1 + 1/5 

2o Thevin equivalent about resistor for piezoceramic only 
%num= [0 1 0]; % 1 over 

oden- [R*L 1 К/С %05+1+ 1/5 


7o Thevin Equivalent for transducer resistor coupling to Silicon layer 

num= [1 0 0] 

den= [R*R2*L*L2 R*L+R2*L2 R*R2*L/C3+R*R2*L2/C1+1 R/C1+R2/C3 
R*R2/(C1*C3)] 


fs= 4400000; 
parper= 1/fs; 
per=1/220000; 
[= |: 


for t=0:parper: 14.2449*per; 
wav(D= (1-exp(-t/(2*per)))*sin(t/per*2*pi); 
ame(D= t: 
I= 1+1; 
end; 


[numd,dend]= bilinear(num,den,fs,220000); 

x= [wav zeros(1,4*length(time))]; 

% Receiving transducer represents a second set of coupling equivalents 
numd2=conv(numd,numd); 


dend2=conv(dend,dend); 


% To see values unconvolved, simply replace numd2 & dend2 w/ numd & dend 


outy= filter(numd2,dend2,x); 
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time2= time+max(time)+1/fs; 
time3= time2+max(time); 
time4= time3+max(time); 
time5= time4+max (time); 


bigtime= [time time2 time3 time4 time5]; 


figure(1) 
w=[100000:100000: 10000000]; 
bode(num,den,w) 


figure(2) 

plot(bigtime,outy, -'); 

title(Ringing for Receiving Transducer Resistor’); 

ylabel(Resistor Response’); 

xlabel( Time (sec)’); 

Zhold on 

Zoplot(bigtime,x,':'); 

Zohold 

To To To To To To To To To To To To To To Fo To To To To To To To To To To To To To To To To To To To To To To To To To To To 
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Transducer Model from HP4195A Freq. Measurements 
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Figure A-2. Direct measurements of transducer frequency response using HP4195A 
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Fig. A-3. Model of response of transducer to input using A-2 
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To Yo Yo Yo To To To To To Yo To To To To Yo To To Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Vo Yo To Ф 
% File saved as hp_ring.m that was used for Figures A-2 and A-3 
clear 
load hp4195a.txt 
i2*pi*hp4195a(:,2); 
MagImp= 1000*hp4195a(:,3); ^ 96Convert kOto Q 
PhaseA= hp4195a(:,4); 
Rpart= MagImp.*cos(pi*PhaseA/180); 
Ipart= MagImp.*sin(pi*PhaseA/180); 
h= Rpart + Ipart.*I; 
nb=4; 
na-5; 
% H is the Complex Freq., f is freq., nb is order of numerator, 
% na is order of denominator, and b and a are real values for 
% the numerator and denominator 


[b,a]= invfreqs(h,f,nb,na); 


w= [1000000: 10000: 10000000]; 
figure(1) 

bode(b,a) 

hold on 

subplot(21 1) 

title(Bode for hp ring.m') 


fs= 4000000; 
parper= 1/fs; 
per=1/220000; 
¡Sl 


for t=0:parper:14.2449*per; 
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wav(I)= (1-exp(-t/per))*sin(t/per*2*p1); 
time(1)= t; 
[= ІЗІ; 


end; 


% Transform descending powers of s to z-transform coefficients 


[bd,ad]= bilinear(b,a,fs); 


x= [wav zeros(1,3*length(time))]; 

% Convert Digital input into Digital Output through the 
% Digital Filter of bd/ad 

outy= filter(bd,ad,x); 


time2= time+max(time)+1/fs; 
time3= time2+max(time ); 
time4= time3+max(time); 


bigtime= [time time2 time3 time4]; 


figure(2) 

bode(b,a,w) 

subplot(211) 

ylabel('! Impedence |") 

title( Transducer Model from HP4195A Freq. Measurements’); 


figure(3) 

plot(bigtime,outy); 

title(Transducer Ringing using HP4195A Freq. Measurements’); 
ylabel( Response to Signal); 

xlabel( Time (sec); 
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figure(4) 
plot(bigtime,x); 

title Transducer Input’); 
ylabel(Signal”); 
xlabel(‘Time (sec)'); 


To To To To To To To To To To To To To Yo To To To To To To To Fo To To To Fo Fo Fo Fo Fo To To To Fo To To To To To To To To 


To Fo To To To To To To To To To To To To To To To To To To To To To To To To To To To To To To To To To To To To To To To To 
%File Saved as Trand2.m that was used for the Bode plots in the Ringing Analysis 


Josection 


clear 

L= 0.980964e-3; 
L2= 0.00084; 

R= 45.5342; 

R2= 35; 

С1= 5.97691е-10; 
С2= 1.23217е-9; 
5-5. 75e-10; 
9: 


W= [1000000:10000:10000000]; 


ишпа=[Г/С2 R/C2 1/(C1*C2)] 
dena=[L R (1/C1+1/C2) 0] 


numb=[L2/C4 R2/C4 1/(C3*C4)] 
denb=[L2 R2 (1/C3+1/C4) 0] 


numS4= L*L2; 
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numS3= L*R2+R*L2; 

numS2= L/C3+R*R2+L2/C1; 

numS1= R/C3+R2/C1; 

numS0= 1/(C1*C3); 

CON= C4+C2; 

denS5= CON*numS4; 

denS4= CON*numS3; 

denS3= CON*numS2+L+L2; 

denS2= CON*numS 1+R+R2; 

denS 1= CON*numS0+1/C14+1/C3; 

denS0= 0; 

bıgNUM= [numS4 numS3 numS2 numS1 numS$0]; 
bigDEN= [denS5 denS4 denS3 denS2 denS1 denSO0]; 


% Plot of the combined parallel equivalent circuit impedance 
% and phase response 

figure(1) 

bode(bigNUM,bigDEN,W) 

hold on 

subplot(211) 

ylabel(Mag. Impedance”) 


% Plot of the impedance and phase of equivalent circuit generated 
% by the HP4195A 

figure(2) 

bode(numa,dena,W) 

hold on 

subplot(211) 

ylabel(Mag. Impedance’) 
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% Plot of the other half of the parallel circuit impedance and phase (alone) 

figure(3) 

bode(numb,denb,W) 

hold on 

subplot(211) 

ylabel( Mag. Impedance’) 

To Fo Fo Yo Yo To To To To To To To Fo Yo To To To To To To To To To Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo 


To To Yo To To To To To To To To Fo Fo To Fo To To To To To To To To To To To To To To To To To To Fo Fo Fo Fo Fo Fo Fo To Fo 
% File saved as ring.m that was used for the input signal and ringing analysis graphs 

% presented in Chapter 3 

clear 

fs= 4000000; 

parper= 1/fs; 

per=1/220000; 

I=1; 


for t=0:parper:14.2449*per; 
wav(I)= (1-exp(-t/per))*sin(t/per*2*pi); 


ите(Г)= t; 
[= 1+1; 
end; 
figure(1) 


plot(time,wav); 


L= 0.980964е-3; 
L2= 0.00084; 

R= 45.5342: 

R2= 35; 

С1= 519769 le-10; 
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(= 1 23217е-9; 

C3= 5.75e-10; 

C4- 2e-9; 

w= [1000000:10000: 10000000; 


numS4= L*L2:; 

numS3= L*R2+R*L2; 

numS2= L/C3+R*R2+L2/C1; 

numS1= R/C3+R2/C1; 

numS0= 1/(C1*C3); 

CON= C4+C2; 

denS5= CON*numS4; 

denS4= CON*numS3; 

denS3= CON*numS2+L+L2; 

denS2= CON*numsS 1+R+R2; 

denS1= CON*numS0+1/C1+1/C3; 

denS0- 0; 

num- [numS4 numS$3 numS2 numS1 numS0]; 
den= [denS3 denS4 denS3 denS2 denS1 4еп50]; 


% Transform descending powers of s into z-transform coefficients 
[numd,dend]- bilinear(num,den,fs); 


x- [wav zeros(1,3*length(time))]; 


% Convert Digital input into Digital Output through the 
% Digital Filter of numd/dend 


outy= filter(numd,dend,x); 


time2= time+max(time)+1/fs; 


time3= time2+max(time); 
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time4= time3+max(time ); 


bigtime= [time time2 time3 time4]; 


figure(2) 
bode(num,den,w) 
subplot(211) 

ylabel(1 Impedance 1); 
title Transducer Model; 


figure(3) 
plot(bigtime,outy); 
title(Transducer Ringing’); 
ylabel('Response to Signal’); 
xlabel("Time (sec)'); 


figure(4) 

plot(bigtime,x); 

titleCTransducer Input‘); 

ylabel( Signal’); 

xlabel(Time (sec)'); 

To To To To To To To To To To To To To To To To Fo To To To To To To Fo To To To To To To To To Yo To To To To To To To To Fo 
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APPENDIX B 


Timing and Burst Generation Board 
B-1 Left half of timing and burst generator schematic 


B-2 Right half of timing and burst generator schematic 


The timing and burst generation board for an unmodified BASS is presented in this 
appendix. Because of the differences between air and water in group velocity and mean 
flow rate, it was necessary to make changes to the initial board. The major modifications 


to the board are detailed in Chapter 3 of the text. 
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Figure B-1. Left half of timing and burst generator schematic 
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APPENDIX C 


Windtunnel Programs 

This appendix contains the programs used with the windtunnel and pitot tube. The 
most interesting result of this section was the difficulty in developing useable libraries in C 
for QBasic. If similar problems are encountered, the shell command of Qbasic coupled to 
a stand alone executable is an excellent solution. Using the programs below, a 286 
computer was capable of taking several measurements per second. Since the response 
time of the pitot tube was 0.5 seconds, the processing speed requirements were exceeded. 
If faster measurements would have been required, a 386, 486, or Pentium could do the 


job. 
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"Qbat.bat - used to load A/D Board libraries so that other libraries could be loaded later 


“inside Obasic. This was not needed in the final version of Windshl.bas 


cd c:\gb45 
ал с1о8 


e сь «өз е» Qua) «ж к= = шь шь = == = өл» өн» өм» «м» MMMM am a am am am am am am am am ае» «ж» eo = 


"Windshl.bas - version of Windtun.bas used to shell to C compiled executables after 


* multiple library development attempts failed 


DECLARE SUB PutBaro (NowBaro AS DOUBLE, AvBaro AS DOUBLE, VarBaro AS 
DOUBLE, SBaro AS DOUBLE, SquareBaro AS DOUBLE, J AS INTEGER) 
DECLARE SUB PutTemp (NowTemp AS DOUBLE, AvTemp AS DOUBLE, VarTemp 
AS DOUBLE, STemp AS DOUBLE, SquareTemp AS DOUBLE, J AS INTEGER) 
DECLARE SUB PutWind (NowSpd AS DOUBLE, Winds) AS DOUBLE, AvSpd AS 
DOUBLE, VarSpd AS DOUBLE, SSpd AS DOUBLE, SquareSpd AS DOUBLE, J AS 
INTEGER) 

DECLARE SUB Temperature (CQ AS INTEGER, NowTemp AS DOUBLE) 
DECLARE SUB Barometer (CQ AS INTEGER, NowBaro AS DOUBLE) 

DECLARE SUB AtoD (DGQ AS DOUBLE, DQ AS INTEGER, COQ AS INTEGER) 
DECLARE SUB Initboard (DQ AS INTEGER) 

DECLARE SUB OutNow (NowBaro AS DOUBLE, AvBaro AS DOUBLE, VarBaro AS 
DOUBLE, NowTemp AS DOUBLE, AvTemp AS DOUBLE, VarTemp AS DOUBLE, 
NowSpd AS DOUBLE, AvSpd AS DOUBLE, VarSpd AS DOUBLE) 

DECLARE SUB SetScreen () 

DECLARE SUB WindSpd (NowBaro AS DOUBLE, NowTemp AS DOUBLE, NowSpd 
AS DOUBLE) 
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‘PROGRAM WINDTUN.BAS 

' This program will include interaction between the subroutines 
‘called Temperature and Barometer and an A/D Board. It will 

' reference a C program that makes use of the MKS Barotron/270 
' Signal Conditioning Package. The A/D board will be use 


"-/- 5 volts for the counts routines. 


! 


First, initialize all variables and arrays that will be used ın your рго- 


t 


gram including the following, which are used in the CALL routines. 


-- a m m m m l a a iu iu au dE E AU SO SO SO SO AUD SB SO AU AUD AUD AU AG cum еш» еқ» әне ан» = = ан» aa ew —— ен» аже аж» ен» ане аж» жне = = = cum Gh a = = че = че a = an Gh m a ED ARD AU ARD A A ӘР A ӘР ҸӘ ӘР ӘР ӘР 


DIM Directions(13) AS INTEGER 'Array for parameters used by 
'СТО8 subroutine CALL.(13 is the 
'maximum number of parameters 
'used in any mode.) 
COMMON SHARED Directions) AS INTEGER ‘Allow subroutine access to 
‘this variable. 
DECLARE SUB CIO8 (MD%, BYVAL DUMMY%, F%) 'Declare subroutine and 
pass 
required parameters. 


F% = 0 ‘Set error flag to 0. 


ON KEY(2) GOSUB Stops 


CONST ADchns = 4 
CONST SIZE = ADchns + 1 
CONST SIZEI = 1200 
CONST VoltsToPU = -9 
CONST bcoef0 = 899.94 
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CONST bcoefl = .03663003# 
CONST rcoef0 = -242.8382 

CONST rcoefl = 2.275128 

CONST rcoef2 = .001879644# 
CONST rcoef3 = -.000004554426# 
CONST rcoef4 = .00000001 132138# 
CONST rcoef5 = -8.142306D-12 


TYPE INFO 
DATESTART AS STRING * 10 
Timestart AS STRING * 8 

END TYPE 


DIM J AS INTEGER 

DIM L AS INTEGER 

DIM COUNTS(1 TO 13) AS INTEGER 
DIM Dataget(1 TO SIZE) AS DOUBLE 
DIM NowTemp AS DOUBLE 

DIM AvTemp AS DOUBLE 

DIM VarTemp AS DOUBLE 

DIM STemp AS DOUBLE 

DIM SquareTemp AS DOUBLE 

DIM NowBaro AS DOUBLE 

DIM SBaro AS DOUBLE 

DIM SquareBaro AS DOUBLE 

DIM AvBaro AS DOUBLE 

DIM VarBaro AS DOUBLE 

DIM Winds(1 TO SIZE1) AS DOUBLE 
DIM NowSpd AS DOUBLE 

DIM SSpd AS DOUBLE 
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DIM SquareSpd AS DOUBLE 
DIM AvSpd AS DOUBLE 
DIM VarSpd AS DOUBLE 
DIM HellfreezesOver AS STRING * 5 
KEY(2) ON 

J=0 

CLS 

HellfreezesOver = "False" 
CALL Initboard(Directions()) 
CALL SetScreen 

SHELL "mainmks.exe" 
Timestart - TIMER 


LOCATE 21,2 
PRINT "Getting first 10 seconds of data: " 


DO 
CALL AtoD(Dataget(), Directions(), COUNTSQ) 
CALL Temperature(COUNTS(), NowTemp) 
CALL Barometer(COUNTS(), NowBaro) 
J=J+1 
CALL PutBaro(NowBaro, AvBaro, VarBaro, SBaro, SquareBaro, J) 
CALL PutTemp(NowTemp, AvTemp, VarTemp, STemp, SquareTemp, J) 


Countup = TIMER - Timestart 
LOCATE 21, 35 
PRINT STR$(10 - Countup) 


LOOP UNTIL Countup >= 10 


[1 





LOCATE 21,2 
ERINT " 


ро 
SHELL "pitotram.exe" 
CALL AtoD(Dataget(), Directions(), COUNTS()) 
CALL Temperature(COUNTSO, NowTemp) 
CALL Barometer(COUNTSO, NowBaro) 
J=J+1 
= 1 
CALL PutBaro(NowBaro, AvBaro, VarBaro, SBaro, SquareBaro, J) 
CALL PutTemp(NowTemp, AvTemp, VarTemp, STemp, SquareTemp, J) 
CALL WindSpd(AvBaro, AvTemp, NowSpd) 
CALL PutWind(NowSpd, Winds(), AvSpd, VarSpd, SSpd, SquareSpd, L) 
CALL OutNow(NowBaro, AvBaro, VarBaro, NowTemp, AvTemp, VarTemp, 
NowSpd, AvSpd, VarSpd) 
LOOP UNTIL HellfreezesOver = "True" 
Stops: 
KEY(2) OFF 
STOP 
END ‘Main 


SUB AtoD (DGQ AS DOUBLE, DQ AS INTEGER, CQ AS INTEGER) 


MD% =4 

FOR 1% = 1 TO ADchns 
CALL CIOS(MD%, VARPTR(D(0)), F%) — ‘Call the subroutine 
DG(I%) = D(O) / 2048 * 5 
С(1%) = 0(0) 

МЕХТ 1% 


ne 








END SUB 'AtoD 


SUB Barometer (CO AS INTEGER, NowBaro AS DOUBLE) 


BCOUNTS = С) + 2048 

IF BCOUNTS < 1 THEN 
BCOUNTS = 1 

ELSEIF BCOUNTS > 4095 THEN 
BCOUNTS = 4095 

END IF 


'The counts values are converted to Barometric Pressure 


NowBaro = bcoef0 + bcoefl * BCOUNTS 


END SUB ‘Barometer 


SUB Initboard (DQ AS INTEGER) 


i Initialize CIO-DASOS using Mode 0 


' Define the variables for MODE (MD%), Base Address (D%(0)) and Range 

' (D%(1)) selected for the CIO-DAS08 and the BUS speed (D%(2)) of your PC. 
' The Base Address is determined by the DIP switch on the CIO-DASOS. 

' The Range is selected by another DIP switch on the CIO-DASOS. 

' The BUS speed can be determined by using the program called PCLK.EXE that 


1 


is included with these example programs. 


m an an an an am an am an an an am еш» ене SO SO SS SOS SCF ss ese «м» өш» «к» «к» ех» «ж» «м» «ше «к» «м» «м» ees «м» ен» «м» «м» е» он» «ж» «м» «м» ен» еш» еш» қи» «м» бы» es cum cum qp өне SB еи» «шә өн ен» өн» UD ене өше «жә өше е» 


MD% = 0 ‘Set MODE to 0 (initialize CIO-DAS08) 
D(0) = &H330 ‘Address as shipped from factory. 


> 





'340 Hex (832 decimal) 

ШІу- 0 'A/D range 0 to 10VDC. 

Do) = 31 ‘Set the time constant for the PCLK 
'31 is typical for an AT type computer 
(30 - ӨМН? bus). 
"This variable is used by Mode 50 
‘and 51, 


CALL CIO8(MD%, VARPTR(D(0)), F%) ‘Call the subroutine 


IF F% <> 0 THEN 
PRINT "Initialization Error” 
STOP 

END IF 


The CIO-DASOS is now initialized. Other modes may now be run. 
' Set upper and lower scan limits of the CIO-DASOS using Mode 1. 

' Define the variables for MODE (MD%), Lower Scan Limit (D%(0)) and Upper 
' Scan Limit (D%(1)) that you wish to select for the CIO-DASOS. This 


example will scan all channels (0 - 1). 


MD% = 1 ‘Set MODE to 1 (Mode to set MUX up) 
D(0) 20 'Lower scan limit 
D(1) = ADchns - 1 ‘Upper scan limit 


‘All channels will be read 
CALL CIO8(MD%, VARPTR(D(0)), F%) — ‘Call the subroutine 
IF F% <> 0 THEN 
PRINT "Channel Initialization Error" 
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STOP 
END IF 


END SUB 'Initboard 


SUB OutNow (NowBaro AS DOUBLE, AvBaro AS DOUBLE, VarBaro AS DOUBLE, 
NowTemp AS DOUBLE, AvTemp AS DOUBLE, VarTemp AS DOUBLE, NowSpd AS 
DOUBLE, AvSpd AS DOUBLE, VarSpd AS DOUBLE) 

OPEN "WindOut.dat" FOR APPEND AS #2 

WRITE #2, TIMER, NowSpd, NowBaro, NowTemp 

CLOSE #2 

LOCATE 6, 65 

PRINT USING "\ \; ТІМЕ$ 

LOCATE 8, 65 

PRINT USING “##.##"; NowSpd 

LOCATE 9, 65 

PRINT USING "##.##"; AvSpd 

LOCATE 10, 65 

PRINT USING "##.##"; VarSpd 

LOCATE 12, 65 

PRINT USING "####"; NowBaro 

LOCATE 13, 65 

PRINT USING "HHHH"; AvBaro 

LOCATE 14, 65 

PRINT USING "####": VarBaro 

LOCATE 16, 65 

PRINT USING "###.##"; NowTemp 

LOCATE 17, 65 

PRINT USING "###.4##", AvTemp 

LOCATE 18, 65 
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PRINT USING "###.##"; VarTemp 


END SUB 'OutNow 


SUB PutBaro (NowBaro AS DOUBLE, AvBaro AS DOUBLE, VarBaro AS DOUBLE, 
SBaro AS DOUBLE, SquareBaro AS DOUBLE, J AS INTEGER) 


SBaro = SBaro + NowBaro 

SquareBaro = SquareBaro + NowBaro * NowBaro 
AvBaro = SBaro / J 

VarBaro = (SquareBaro / J) - AvBaro * AvBaro 


END SUB 'Ри Ваго 


SUB PutTemp (NowTemp AS DOUBLE, AvTemp AS DOUBLE, VarTemp AS 
DOUBLE, STemp AS DOUBLE, SquareTemp AS DOUBLE, J AS INTEGER) 


STemp = STemp + NowTemp 

SquareTemp = SquareTemp + NowTemp * NowTemp 
AvTemp = STemp / J 

VarTemp = (SquareTemp / 7) - АУТетр * АУТетр 


END SUB 'РиТетр 


SUB PutWind (NowSpd AS DOUBLE, Winds() AS DOUBLE, AvSpd AS DOUBLE, 
VarSpd AS DOUBLE, SSpd AS DOUBLE, SquareSpd AS DOUBLE, J AS INTEGER) 


' Winds(J) = NowSpd 
SSpd =SSpd + NowSpd 
SquareSpd = SquareSpd + NowSpd * NowSpd 
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AvSpd = SSpd /J 
VarSpd = (SquareSpd / J) - AvSpd * AvSpd 


END SUB 'PutWind 


SUB SetScreen 
LOCATE 4, 51 
PRINT "DATE: " 
LOCATE 4, 57 
PRINT USING \ \; DATES 
LOCATE 5, 51 
PRINT USING " Start:\ \"; TIME$ 
LOCATE 6, 51 
PRINT "Current Time: " 
LOCATE 8, 51 
PRINT "Wind Speed: m/s" 
LOCATE 9, 51 
PRINT "Running Mean: m/s" 
LOCATE 10, 51 
PRINT "Run Variance: m/s" 
LOCATE 12, 51 
PRINT "Pressure: mBars" 
LOCATE 13, 51 
PRINT "Running Mean: mBars" 
LOCATE 14, 51 
PRINT "Run Variance: mBars" 
LOCATE 16, 51 


PRINT "Temperature: С 
LOCATE 17,51 
PRINT "Running Mean: С" 
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LOCATE 18.51 
PRINT "Run Variance: С" 


END SUB 'SetScreen 


SUB Temperature (CQ) AS INTEGER, NowTemp AS DOUBLE) 


Yoh 2 2 2 K k k k k K k kk kK DK ck ok K New Temperature Probe ЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖЖ 


' This Temperature Probe usually gives values between 2.7V and 3.1V to 
‘the A/D Board. 2.7V is 0 degrees C. 3.1V is 


NowTemp = C(2) / 2048 * 5 * 100.4 - 273.16 


END SUB ‘Temperature 


SUB WindSpd (NowBaro AS DOUBLE, NowTemp AS DOUBLE, NowSpd AS 
DOUBLE) 


DIM AirDens AS DOUBLE 
DIM Diff AS DOUBLE 


OPEN "DNDiff.tmp" FOR INPUT AS #1 
INPUT #1, Diff 
CLOSE #1 


AirDens = .34838 * NowBaro / (NowTemp + 273.15) 
NowSpd = 16.306 * (SOR(Diff / AirDens)) 


END SUB 'WindSpeed 
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/* MainMKS.C - version of MKSSet.CPP used to create the executable MainMKS.exe*/ 


/* which is used in the program Windshl.bas*/ 


/* Filename: MKSSet.CPP */ 
#include <graph.h> 
#include <stdio.h> 
#include <stdlib.h> 
#include <conio.h> 
#include <dos.h> 
#include <bios.h> 
#include <string.h> 


#include <malloc.h> 


/* PARALLEL INTERFACE FOR METRABYTE (8255) */ 
#define PORTA 0x300 
#de fine PORTB PORTA + 1 
#define PORTC PORTA + 2 
#define PPI CTRL PORTA + 3 


/* DECLARATION OF SUBROUTINES */ 


void hard_init(void); 


main () 
| 
/* Set up the Parallel Metrabyte Hardware */ 
hard init(); 


no 





| /* End of Main Program  */ 


void hard_init(void) 
( 
/* Set the PPI board for Port A, Port B, Port C-lower (INPUT) 
and Port C-upper (OUTPUT). Control word #11= 0x93 */ 
outp(PPI_CTRL,0x93); 


/* Release the HOLD\ line with Port C bit set */ . 
outp(PPI_CTRL,0x09); 


/* Test DECIMAL OP\line to insure MKS display unit is set right */ 
if (inp(PORTC) & 0x04) 
{ 
/* Line is High (Decimal Point is WRONG) - warning */ 
printf(^nDisplay decimal point out of position - check MKS 270C setup. \n"); 


/* Waits for Keystroke and then Clear Screen to Continue */ 
printf(^nHit Any Key When Ready\n"); 

getchQ; 

_clearscreen(_GCLEARSCREEN); 

} 


} /* End of Function hard_init */ 


/* Pitotram.C - version of Pitot.cpp used to create the executable Pitotram.exe*/ 
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/* which ıs used in the program Windshl.bas*/ 


/* Program Pitot.cpp */ 
#include <graph.h> 
#include <stdio.h> 
#include <stdlib.h> 
#include <conio.h> 
#include <dos.h> 
#include <bios.h> 
#include <string.h> 


#include <malloc.h> 


/* PARALLEL INTERFACE FOR METRABYTE (8255) */ 
#define PORTA 0x300 
#define PORTB PORTA + 1 
#define PORTC PORTA + 2 
Hdefine PPI_CTRL PORTA + 3 


/* DECLARATION OF SUBROUTINES */ 
void pitot_diff(void); 


/* Globals for the Program */ 
double Diff; /* Primary data value */ 
PILCE *fp; 


void main!) 


{ 
/* Get Diff using Parallel Metrabyte Hardware */ 
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pitot. diff(); 


| /* Endof Main Program  */ 


void pitot_diff(void) 
{ 
unsigned char msb,hi о; 
fp= fopen("DADIff. tmp”, “w”); 


/* Activate the HOLD\ line with Port C bit clear */ 
outp(PPI_CTRL,0x08); 


/* Get the BCD reading from the Pitot Tube */ 
msb = inp(PORTC) & 1; 
hi = inp(PORTB), 
lo = inp(PORTA); 


/* Release the HOLD\ line with Port C bit set */ 
outp(PPI_CTRL,0x09); 


Diff = (double)msb + (double)hi/100 + (double)lo/10000; 
fprintf(fp,"%lf\n" Diff); 
fclose(fp); 

} /* End of Function pitot_diff */ 
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APPENDIX D 


TTBasic Program for BASS 


The TTBasic program which follows has been modified from ıts usual format. The 


delay between measurements has been increased to two seconds by the command ‘sleep 


200", and the program is set to read only one path of data. 
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10 PRINT " Single Axis Air 7-18-95 from BigVRamC.TTB RAM" 
20 Q=&H4000 :REM DATAFILE STARTS AT 4000 IN RAM 
20 Q=&H0143 :REM DATAFILE STARTS AT 0143 IN ROM 


100 X=&H112 

130 ASM X,JSR & HFFE2 :REM FLUSH UART BUFFER 

140 ASM X,CLR &H11 = :REM DISABLE UART 

141 ASM X,LDAA #5;STAA &H10:REM BAUD RATE 

142 ASM X,OIM &H42,&H1B :REM ENABLE INTS 

143 ASM X,LDAA #95;STAA &H1C:REM GET 100 Hz RATE 

144 ASM X,OIM &H1B,&H11 :REM RESTART UART 

150 ASM X,RTS 

160 CALL &H112,0 -REM FLUSH UART BUFFER 

200 ASM &H9E,DB &H00 :REM CONTROL DISABLE 

202 SLEEP 0:PCLR 0,1,2,3,4,5,6,7,8,14,15:PSET 16:REM ESTABLISH DDRs 
208 SLEEP 240:PCLR 16:REM LET A/D CALIBRATE FINISH IF STARTED. 
210 PRINT" TYPE WAKE TO GET CONTROL (5sec)" 

212 X = 0: STORE X,#4,&H01010101 : X 2 0: ITEXT X,500 

214 X = 0: IF GET(X,#4)<>&H57414B45 GOTO 300 


230 STOP 

300 ONERR 100 

‘ASSEMBLY ROUTINES 

1010 X=0:A=0:B=0:C=0:D=0:E=0:F=0:G=0:H=0:K=0:L=0:M=0:N=0 
1015 PRINT #6H,A,B,C,D,E,F,G,H,K,L,M,X 

1020 GOSUB 8000:REM FIRST PASS 

1025 PRINT #6H,A,B,C,D,E,F,G,H,K,L,M,X 

1030 GOSUB 8000:REM SECOND PASS 

1031 PRINT #6H,A,B,C,D,E,F,G,H,K,L,M,X 
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' CALIBRATE A/D 

PODOS CEEP 0:PSET 15:SLEEP 600:PSET 16:REM PWR ON 

2070 PCLR 15:SLEEP 240:REM A/D CALIBRATE;2,882,040 CY CLES АТ 1.2288 
MHz 

2090 N=0 


3000 RTIME 

3040 X=0:REM INITIALIZE DATAFILE 

3050 STORE X,#1,&HEE:REM ESTABLISH STRAD AT ORIGIN 

3055 STORE X,#1,?(2):REM HOURS 

3060 STORE X,#1,?(1):REM MINUTES 

3070 STORE X,#1,?(0):REM SECONDS 

3075 STORE X,#1,N:REM COUNT 

3090 SLEEP 0 

3100 С-“%Н73000000-О--Х:КЕМ MUX ADDRESS AND DATAFILE POINTER 
‘AB IS ADDRESS OF FLAGGED CONVERSIONS 

3110 CALL A,C,X :REM CALL A/D ROUTINE,POWER LEFT ON AT END 
3111 X=X%&H10000-Q:REM GET BACK THE DATAFILE POINTER 
‘SUBTRACT ROUTINE (CALL B) EXPECTS FLAGGED DATA AT 6F38. 
3119 C=&H73000000+Q+X:REM DATAFILE LOCATION 

3120 CALL B,C,X:REM CALL SUBTRACT AND TRANSFER 

3121 X=X%&H10000-Q:REM GET BACK THE BASIC DATAFILE POINTER 
3195 SLEEP 200 

3200 N=N+1 


3210 X=0 

3220 PRINT #2H,GET(X #1) 

3230 PRINT 402,GET(X,1),":" ,GET(X,&1),":" ,GET(X,£1)," :REM HR:MM:SS 
3240 PRINT #02, GET(X,#1) 

3260 PRINT #04H,GET(X,#2) 
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3270 PRINT 
'3280 PRINT {0,6};:REM OUTPUT AS BINARY FOR SPEED 
3300 GOTO 3000 


8000 X=&H7300:REM MULTIPLEXOR LIST 

8010 ASM X,DW £H00FF;¡DW £H01FF;¡;DW &HFFFF;DW &HFFFF:REM ACM 1/A 
‘8015 ASM X,DW &H0888;DW &H0989;DW &H0A8A;DW &HOB8B:REM ACM 
1/A5-A8 

'8020 ASM X,DW &H20FF;DW &H21FF;DW &H22FF;DW &H23FF:REM ACM 2 
'8025 ASM X,DW &H28FF;DW &H29FF;DW &H2AFF;DW £H2BFF 

'8030 ASM X,DW &H4000;DW &H4100;DW &H4200;DW &H4300:REM ACM 3 
'8035 ASM X,DW &H4800;DW &H4900;DW &H4A00;DW &H4B00 

'8040 ASM X,DW &H6000;DW &H6100;DW &H6200;DW &H6300:REM ACM 4 
'8045 ASM X,DW &H6800;DW &H6900;DW &H6A00;DW &H6B00 

:8050 ASM X,DW &H8000;DW &H8100;DW &H8200;DW &H8300:REM ACM 5 
‘8060 ASM X,DW &H8800;DW &H8900;DW &H8A00;DW &Н8В00 

'8070 ASM X,DW &HA000;DW &HA100;DW &HA200;DW &HA300:REM ACM 6 
‘8080 ASM X,DW &HA800;DW &HA900;DW &HAA00;DW &HABOO 

8082 ASM X,DW &HC000;DW &HC100;DW &HC200;DW &HC300:REM ACM 7 
$084 ASM X,DW &HC800;DW &HC900;DW &HCA00;DW &HCB00 

'8086 ASM X,DW &HE000;DW &HE100;DW &HE200;DW &HE300:REM ACM 8 
8088 ASM X,DW &HE800;DW &HE900;DW &HEA00;DW &HEBOO 

8090 ASM X,DW &HFF00;DW &HFF00:REM END OF LIST 


8100 X=&H3800:A=X:REM A/D ROUTINE 

"X IS POINTER TO AUX ARRAY, AB IS POINTER TO FLAGGED CONVERSION 
STORAGE 

8105 ASM X,STD &H7000:REM POINTER TO DATAFILE 

8107 ASM X,LDD #&H6F38;STD &H7002:REM POINTER TO FLAGGED 
CONVERSIONS 
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8110 ASM X,OIM &H04,&H17;PSHX;LDX #&H0600:REM POWER ON, WAIT 5ms 
8120 DZX:ASM X,DEX;BNE D;PULX:REM WAIT (4*.81us/LOOP) 

8128 ASM X,SEI:REM DISABLE INTERRUPTS FOR TIMING 

8150 ASM X,CLR &H11;LDAA #&H10;STAA &H10:REM DISABLE UART, SET 
SCI FOR 

8160 ASM X,OIM 8,&H11:REM INTERNALLY CLOCKED 8 BIT DATA RECEIVE 
ENABLE 

‘DUMMY PULSE SENT OUT TO INITIALIZE 

8170 ASM X,OIM &HAB,&H15;AIM &HBF,&H15:REM DUMMY AUX LOCATION 
8171 ASM X,OIM &H80,&H03;0IM &H80,&H03:REM START TIMING P27=1 

8172 ASM X,AIM &H7F,&H03:REM REMOVE PULSE P27=0 

8174 D=X:ASM X,LDAB #&H02;BITB &H03;BEQ D:REM CHECK P21 FOR DONE 
8179 C=X:ASM X,AIM &H14,&H15:REM ENTRY POINT FOR AUX,MASK PORT 5 
8180 ASM X,LDAA £H00,X;INX;INX:REM TOP OF MUX LIST,MOVE TO NEXT 
8185 ASM X,ORAA &HI5;STAA &H15:REM PUT MUX ON PORT 5 

8200 ASM X,OIM &H80,&H03;0IM &H80,&H03:REM START TIMING P27=1 

8210 ASM X,AIM &H7F,&H03:REM REMOVE PULSE P27=0 


8220 D=X:ASM X,LDAB #&H02:REM TEST FOR DONE 
8230 ASM X,BITB &H03:REM CHECK P21 
8240 ASM X,BEQ D 


8250 ASM X,LDAA &H00,X:REM GET NEXT MUX WORD 

8260 ASM X,STAA &H15:REM PUT MUX ON PORT 5 

8270 ASM X,INX;INX:REM MOVE MUX POINTER TO NEXT WORD 

8280 ASM X,OIM &H80,&H03;0IM &H80,&H03:REM START TIMING P27=1 
8290 ASM X,AIM &H7F,&H03:REM REMOVE PULSE Р27=0 


8400 ASM X, CPX 0; CPX 0; CPX 0; CPX 0 
8402 ASM X, CPX 0; CPX 0; CPX 0; CPX 0 
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8403 АЗМ Х, СРХ 0; СРХ 0; СРХ 0; СРХ 0 
8404 АЗМ Х, СРХ 0; СРХ 0; СРХ 0; СРХ 0 


8410 ASM X,LDAA &H11; LDAB &H12:REM TRIGGER A READ 

8420 ASM X, CPX 0; CPX 0; CPX 0; CPX 0:REM 16ECLOCK 

8430 ASM X, LDAA &H11;LDAB &H12:REM HIGH BYTE IN B AND READ NEXT 
8440 ASM X, CPX 0; CPX 0; CPX 0; CPX 0:REM 16ECLOCK 

8450 ASM X, LDAA &H12:REM LOW BYTE IN A DON'T READ NEXT 

‘REM FLIP THE BITS OF THE WORDS 

8451 ASM X, ROLA;RORB;ROLA;RORB;ROLA;RORB;ROLA;RORB 

8452 ASM X, ROLA;RORB;ROLA;RORB;ROLA;RORB;ROLA;RORB;ROLA 


8480 ASM X,TIM &H20,&H03:REM TEST BOTH RECEIVED 

8485 ASM X,BEQ G:REM BRANCH IF RECEIVED 

8490 ASM X,CLRA;CLRB:REM IF NOT RECEIVED,CLEAR 

8505 G=X 

‘PUT AWAY TO @(50) ARRAY THE FLAGGED VALUES 

8510 ASM X, PSHX; LDX &H7002; STD 0,X; INX; INX; STX &H7002; PULX 


8520 ASM X,LDAB #&HFF 

8530 ASM X,EORB &H00,X:REM CHECK FOR END OF LIST 
'8540 ASM X,BNE D:REM LOOP IF NOT DONE 

8540 ASM X, BEQ L;JMP D : L=X 

8550 ASM X,LDAB #&H02:REM TEST FOR HOLD PULSE 
8555 D=X 

8560 ASM X,BITB &H03:REM CHECK P21 

8570 ASM X,BEQ D 


8572 ASM X, CPX 0; CPX 0; CPX 0; CPX 0 


88 





8573 ASM X, CPX 0; CPX 0; CPX 0; CPX 0 

8574 ASM X, CPX 0; CPX 0; CPX 0; CPX 0 

8575 ASM X, CPX 0; CPX 0; CPX 0; CPX 0 

8580 ASM X,LDAA &H11; LDAB &H12:REM TRIGGER A READ 

8590 ASM X, CPX 0; CPX 0; CPX 0; CPX 0:REM 16ECLOCK 

8600 ASM X, LDAA &H11;LDAB &H12:REM HIGH BYTE IN B AND READ NEXT 
8610 ASM X, CPX 0; CPX 0; CPX 0; CPX 0:REM 16ECLOCK 

8620 ASM X, LDAA &H12:REM LOW BYTE IN A DON'T READ NEXT 

‘REM FLIP THE BITS OF THE WORDS PUTS HIGH BYTE IN A, LOW BYTE IN B 
8621 ASM X, ROLA;RORB;ROLA;RORB;ROLA;RORB;ROLA;RORB 

8622 ASM X, ROLA;RORB;ROLA;RORB;ROLA;RORB;ROLA;RORB;ROLA 


8630 ASM X,TIM &H20,&H03:REM TEST BOTH RECEIVED 
8635 ASM X,BEQ H:REM BRANCH IF RECEIVED 

8645 ASM X, CLRA;CLRB:REM IF NOT RECEIVED,CLEAR 
8655 Н=Х 


‘FINISH PUT AWAY TO @(50) ARRAY THE FLAGGED VALUES 
8680 ASM X, PSHX; LDX &H7002; STD 0,X; INX; INX; STX &H7002; PULX 


8705 ASM X,CLR &H11:REM DISABLE SCI 

8710 ASM X,LDAA #5;STAA &H10:REM BAUD RATE FROM TIMER] (#5 GIVES 
9600) 

8720 ASM X,OIM &H42,&H1B:REM ENABLE TIMER2 INTS AND E/128 

8730 ASM X,LDAA #95;STAA &H1C:REM TIMER COUNT TO GET 100 Hz 

8740 ASM X,OIM &H02,&H11:REM 02 RESTART ПАКТ NO RECEIVE 

‘reconf uart, ENABLE IT, power off (PHYSICALLY ENABLE IT) 

8748 ASM X,CLI:REM CLEAR INTERRUPT MASK 
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8750 ASM X,AIM &H14,&H15:REM PUT MULTIPLEXORS ON PARK 
8760 ASM X,AIM &HFA,&H17:REM P60,P62=0,AUX & POWER OFF 
8775 ASM X,LDD &H7000:REM NO AUX VALUES STORED 

8780 ASM X,RTS 


9100 B=X:REM SUBTRACT AND TRANSFER SUBROUTINE USING STRMEM 
‘POINTERS: m(7000) = DATA ARRAY POINTER (NOT UPDATED WITH 
STRMEM) 

' m(7002) = FLAGGED AXIS PAIRS STORAGE 

' 6F38 IS @(50) INRAM MUST BE SAME AS A/D ROUTINE 

' m( X )= AUX LIST FOR END OF LIST CHECKING. PASSED IN X 


9110 ASM X, STD &H7000 

9115 ASM X, LDD #&H6F38; STD &H7002 

9120 K=X: ASM X,PSHX;LDX &H7002;LDD 0,X;INX;INX:REM GET WORD 
9140 ASM X,BEQ F:REM TEST FLAG ON NORMAL MEAS 


9150 ASM X,SUBD 0,X:REM DOUBLE SUBTRACT 
9152 ASM X,RORA;RORB:REM SHIFT RIGHT WITH CARRY = DIVIDE BY TWO 


9160 ASM X,TST 0,X:REM TEST FLAG ON REVERSED MEAS 
9170 ASM X,BEQ F 


9200 DZX:ASM X,INX;INX;STX &H7002;PULX:REM MOVE ON TO NEXT PAIR 


9203 REM 'SAVE RESULT TO DATAFILE WITH STRMEM 


9205 ASM X,PSHX;LDX &H7000;STD O,X;INX;INX;STX &H7000;PULX 
9205 ASM X,PSHB;JSR &HFFD3;PULA;JSR &HFFD3 
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' CHECK FOR END OF LIST. AUX LIST= FOUR BYTES / AXIS 
9240 ASM X,INX;INX;INX;INX 

9250 ASM X,LDAB #&HFF:REM CHECK FOR END OF LIST 
9255 ASM X,EORB &H00,X 

9260 ASM X,BNE K:REM LOOP IF NOT DONE 


0265 ASM X,LDD &H7000:REM LOAD DATAFILE POINTER FOR RETURN 


9270 ASM X,RTS:REM EXIT 


9275 F=X 


REM FLAG MISSED 


9280 ASM X,LDD #&H8000 


9290 ASM X,BRA D 


9900 RETURN 


9999 END 


> 
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APPENDIX E 


Specifications of Manufacturer for E-188/220 


This appendix contains the manufacture’s specifications for the E-188/220 transducer 
produced by Massa Products. The agreement that was produced between these 
specifications and the results from the HP4195A spectrum analyzer was outstanding. The 
measurements of four different E-188/220 transducers using the HP4195A produced good 


agreement between the transducers and the specifications. 
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: Model E:188 


Specifications 


`..” Е-188/215 


Frequency at. | | 
Receiving 
Sensitivity .. 
Bandwidth- 
(Transmitting) 
Transmitting 
Sensitivity 

(dBvs 1ubar) 
pervolt at. 

1 foot) 
Receiving : 
Sensitivity . 
(dB vs 1 voit/ubar) 


Driving Voltage | 
(10% duty cycle) ~. 


Nominal Ітредепсе 


Tota! Beam Angie | 
(-3dB) Conical 


Operating . 


Temperature ^ ^ — 


Humidity | 


E-188/220 


220kHz 
z25kHz . 


215kHz - 
z2.5 kHz 


25kHz 


25kHz e 


+20 


sov Р-Р 
100092 10000 
Т 10°. 


0-70°С 0-70°С 


0-90% Non-oondensing D 


E = 1VRMS 


98 $ 1ыБаг ат М 


OUTLINE DIMENSIONS 





















580 2400 - 50 
$10 DIA — = 
114 n (009 6) 
50 
E ‚ 0! ¡e 
13 (12 7) 
480 DIA (3 3) СОО. DIA NOM 25 
(12 191 | 16 41 
| (2 5! 









RTI 


/ eom 
/ ALUM HOUSING ( MILLIMETERSI 
ACOUSTIC FACE FREQUENCY IDENTIFIER 


WHITE DOT 21ShHe > 2kg 
BROWN DOT 220kHz * 2xHrz 


TRANSMITTING RESPONSE 





RECEIVING RESPONSE 
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Weight- 10grams . 10grams `: 8 80 
io 
І P р , С 100 
All specifications typical at 22°C and Barometric pres-- 209 a 220 230 240 250 
sure of 751mm of Mercury and subject to change wıth- Frequency «на 
out notice. . a eee 
— IMPEDANCE 





DIRECTIONAL RESPONSE 


is. 
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Massa. Products: 
Corporation 


ATT 
Hingham, Massachusetts 02043 


25. 


Tel:617:749-4800 
TWX: 710-348-6932. 








Figure E-1. Massa Products” specifications for E-188/220 
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